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A model-independent, three-dimensional source function for pion pairs has been 
extracted from Pb+Pb collisions at \/s NN = 17.3 AGeV. The extracted source ex- 
hibits long-range non-Gaussian tails in the directions of the pion-pair net transverse- 
momentum and of the beam. Comparison with the Therminator model allows for an 
extraction of the pion source proper breakup time and of emission duration in the 
collisions. 



1. INTRODUCTION 

A deconfined phase of nuclear matter is expected to be formed at the high energy densities 
created in relativistic heavy ion collisions [1]. It is expected that the signatures of such a 
phase are reflected in the space-time extent and shape of particle emission source-functions. 

Recently, one-dimensional (ID) source imaging techniques [2-4] have revealed a non- 
trivial long-range structure in the two-pion emission source at RHIC [5, 6]. The origins of 
this structure are still investigated. The nature of such structures and the potentially useful 
information they can provide could be revealed, on one hand, by examining the possible 
presence of such a structure in heavy-ion collisions at intermediate SPS energies and, on 
the other, by examining the dependence of a structure on the direction within the source. 
The NA49 Collaboration has investigated Pb+Pb collisions over a wide range of bombarding 
energies at the CERN SPS during the last decade [7], accumulating pion data. The rich data 
set provides a unique opportunity to search for non-trivial structures at the SPS, investigate 
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their dependence on energy and direction and extract any useful information those structures 
may provide. 

In this paper, the first three-dimensional (3D) emission source image for low px pion 
pairs produced in central Pb+Pb collisions at a/%7v = 17.3 GeV is presented. The result 
is compared with calculations from the Therminator model in order to extract the source 
proper breakup time and emission duration. 

2. SETUP AND ANALYSIS 

2.1. Experimental Setup 

The data presented here were taken by the NA49 Collaboration during the years 1996- 
2002. The incident beams of 158 AGeV were provided by the CERN SPS accelerator. The 
NA49 Large Acceptance Hadron Detector [8] achieves a large-acceptance precision-tracking 
(Sp/p 2 ~ (0.3 — 7). 10~ 4 (GeV/c) -1 ) and particle identification using time projection chambers 
(TPC's). Charged particles are detected by the tracks left in the TPC's and identified by the 
energy deposited in the TPC gas. Mid-rapidity particle identification is further enhanced 
by a time-of-flight wall (resolution of 60 ps). Event centrality is determined by a forward 
calorimeter which measures the energy of spectator matter. 

2.2. Three-Dimensional Correlation Function 

The 3D correlation function, C(q), was calculated as the ratio of pair to uncorrelated 
reference distributions in relative momentum q for 7r + 7r + and 7r~7r~ pairs. Here, q = ( pi ~ P2 ) 
is half of the momentum difference between the two particles in the pair Center-of-Mass 
System (PCMS) frame. The pair distribution was obtained using pairs of particles from the 
same event and the uncorrelated distribution was obtained by pairing particles from different 
events. The Lorentz transformation of q from the laboratory frame to the PCMS was done 
by a Lorentz transformation to the pair Locally Co-Moving System (LCMS) frame along the 
beam direction followed by a Lorentz transformation to the PCMS along the direction of the 
transverse momentum of the pair. 

Track merging and splitting effects were suppressed by appropriate cuts on both the pair 
and uncorrelated distributions. The pair cuts require the two particles in each pair to be 
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separated by at least 2.2 cm over 50 pad rows in the vertex TPC's. A 20% increase in this 
minimum separation results in the correlation data points fluctuating within the statistical 
errors. Hence, the systematic uncertainty associated with the pair cuts is deemed smaller 
than the statistical uncertainty. 

The effects of track momentum resolution were assessed by jittering the momentum of the 
tracks in the data by the maximum momentum resolution, 5p/p 2 ~ 7.10~ 4 (GeV/c) -1 , and 
by re-computing the 3D correlation function. The resulting correlation function incorporates 
then twice the effect of the momentum resolution and was found to be very close to the raw 
un-smeared correlation function. This is not surprising, considering that the mean momen- 
tum of the tracks used in this analysis is 1.4 GeV/c, resulting in a momentum resolution of 
5p/p ~ 0.1%. The above procedure gives a very small smear to the correlation peak in the 
region q < 10 MeV. In practice, this tends to reduce the imaged source intensity at large 
separation. As a consequence, the results for the source function presented here may slightly 
underestimate the actual source function at large r. 



2.3. Angular Moments 

In the Cartesian surface-spherical harmonic decomposition technique [5, 9], the 3D cor- 
relation function is expressed as 

C(q) - 1 = R(q) = J2H <•..«« (?)<•..*, (1) 

/ Ot\...Otl 

where / = 0, 1,2, . . ., cej = x,y or z, A l ai ^ (0,^) are Cartesian harmonic basis elements (J7 q 
is solid angle in q space) and R l ai „ mai (q) are Cartesian correlation moments given by 

R l «,..«M) = J ^<...<,(^(q) ( 2 ) 

where q is the modulus of the 4-vector q. The coordinate axes are oriented so that z is 
parallel to the beam (long) direction, x points in the direction of the total momentum of the 
pair in the locally co-moving system (LCMS) (out), and y is perpendicular to the first two 
axes (side). 

The correlation moments, for each order /, may be calculated from the measured 3D 
correlation function using Eq. (2), but can be then vulnerable to directional inefficiencies. 
Alternatively, the moments may be fitted using Eq. (1), avoiding regions of poor efficiency 
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or accounting for those in the fit weights. In following the latter type of analysis, Eq. (1) has 
been truncated at / = 4 and expressed in terms of independent moments only. The higher 
order moments have been found to be negligible. Up to order 4, there are 6 independent 
moments: R°, R\ 2 -, -Ry2? ^x4> RyA an d ^x2j/2 where R\ 2 is shorthand for R\ x etc. These 
independent moments were extracted as a function of q, by fitting the truncated series to 
the measured 3D correlation function with the moments as the parameters of the fit. 



Each independent correlation moment can be imaged using the ID source imaging code of 
Brown and Danielewicz [2-4] to obtain the corresponding source moment within each order /. 
Both the effects of Bose-Einstein symmetrisation and of Coulomb interaction, the sources of 
the observed correlations, are accounted for in the source imaging code. Thereafter, the total 
source function, S(r), can be reconstructed by combining the source moments for each I: 



Another method of construction of the 3D source function from the 3D correlation function 
is by fitting the latter directly upon assuming the shape for the source function. Since the 
3D correlation function is represented by the Cartesian moments in the Cartesian harmonic 
decomposition, this amounts to fitting the six independent moments, dependent on q, with 
a trial source function. Two trial functional shapes were considered in this analysis: the 3D 
Gaussian, frequently termed an ellipsoid shape, and the Hump function, with six adjustable 
parameters, given by 



2.4. Source Reconstruction 




(3) 



S(x,y,z) = Aexp[-/ S (— ) - /,( 




)] 



(4) 



where f s = 1/[1 + (r/r ) 2 )] and /, = 1 - f s . 



3. RESULTS 



Figure 1 (a) shows a comparison between the ID correlation function C(q) and the / = 
moment 1+R(q) for mid-rapidity low-pr pion-pairs from central Pb+Pb collisions. The func- 
tions are in very good agreement with each other as expected in the absence of significant 
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detection inefficiencies with dependence on angle in q-space. The figure also shows that the 
Hump function (solid line) fits the data very well while the ellipsoid fit (dotted line) under- 
estimates the data at q < 13 MeV/c. This difference in correlation fits at low q corresponds 
to a difference in the source function at large r > 15 fm, as evident in Fig. 1 (b). The Hump 
shape (solid line) is in good agreement with the source image (squares) whereas the ellipsoid 
shape (dotted line) underestimates the image. Given that the discrepancy occurs at large r, 
that discrepancy becomes even more pronounced for the radial density in Fig. 1(c). 

Figure 2 shows the different / = 2 and / = 4 moments (open circles) as a function of pion 
separation. The / = 4 moments are significantly smaller in magnitude compared to the / = 2 
moments: this justifies the truncation of the series Eq. (1) at / = 4. The Hump fit (solid 
line), with x 2 /ndf = 1.3, is in close agreement with the data whereas the ellipsoid function 
(dotted line) gives a poor fit to the data (x 2 /ndf = 6.8), as is visually evident. 

Figure 3 (a)-(c) shows the 3D correlation function profiles in the x, y and z directions 
while the corresponding source function profiles are shown in panels (d)-(f). The source image 
(squares) is in good agreement with the Hump function fit (solid line), but is underestimated 
by the ellipsoid fit (dotted line) in all 3 directions. These non-Gaussian tails in the Hump 
source function profiles have corresponding manifestations in the correlation function profiles 
at low q, relative to the ellipsoid fit. 

4. THERMINATOR MODEL 
4.1. Model Assumptions 

The Therminator model [10] can shed more light on the source breakup and emission 
dynamics. The Therminator model incorporates the following: (1) the Bjorken assumption 
of longitudinal boost invariance, (2) blast-wave expansion in the transverse direction with 
transverse velocity profile semilinear in transverse radius p [11], (3) thermal emission of 
particles from a cylinder of infinite longitudinal size and finite transverse radius p max , and 
(4) all known resonance decays. 

Assumptions (1) and (2) of the model imply that the expanding source consists of fluid 
elements shaped like onion rings with their axes aligned with the beam axis. Each fluid 
element expands transversely with a transverse velocity semi-linear in p, as well as translates 
longitudinally with a velocity profile linear in z. 
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Each fluid element breaks up after a proper breakup time r in its own rest frame, with 
subsequent particle emission. Bjorken assumption of longitudinal boost invariance implies 
that the breakup (emission) time, t, of a source element at position z in the laboratory frame 
is given by t 2 = r 2 + z 2 . 

At the point of source breakup, particles within the Therminator model leave each source 
element defined by their longitudinal and transverse positions: z and p. All particle emissions 
are externally viewed as occuring from a freeze-out hypersurface defined in p-r plane as 
r = To + ap, where a is a parameter which controls the space-time correlation and tq is the 
proper breakup time of the source element located at p = 0, i.e. on the beam axis, and r 
is proper breakup time of the source element located at finite p. Hence, particles which are 
emitted from a generic fluid element with coordinates (z,p) will have emission time t in the 
laboratory frame given by t 2 = (r + ap) 2 + z 2 . 

In this analysis, Therminator is used in the blast-wave mode and a is set to the negative 
value of —0.5 [12]. Hence, r decreases linearly with p down to a minimum value at p = p max - 
The negative value of a implies a negative space-time correlation, i.e. particles at large p 
freeze-out earlier. Hence, the source emits or burns from outside in. 

Since each source element is defined by only one value of the proper breakup time r, all 
particle emissions from this fluid element happen instantaneously in the rest frame of the 
source element and the proper emission duration in the above scenario is 0. On the other 
hand, one can allow for a finite non-zero proper emission duration in each source element, 
if the single proper breakup time r of each source element is replaced with a distribution of 
breakup times. This corresponds physically to the source element breaking up over a finite 
time interval rather than instantaneously. 

One such parametrization is an exponential distribution of breakup times r' with a width 
At given by dN/dr' = @ ^ T A ~ T ^ exp[— (V — r)/Ar] . The minimum value of r' is the initial 
single breakup time r for that source element. With such a distribution of breakup times, 
each source element emits particles from a family of hypersurfaces, each of which is defined 
by a t' value which is sampled according to the exponential distribution. In such a scenario, 
particle emission from each source element takes place over a finite time duration defined by 
At rather than instantaneously. In this parametrization, At represents the proper emission 
duration in the rest frame of the fluid element. 
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4.2. Comparison to Data 

Figure 4 shows a comparison of the source image (squares) with the calculated source 
function (circles) from the Therminator model in its blast-wave mode. A good match is 
obtained with values of r = 7.3 fm/c and At = 3.7 fm/ c, when all resonance decays are 
turned on. The good agreement between extracted source image and calculated source func- 
tion indicates that the data are consistent with the basic ingredients of Bjorken longitudinal 
expansion and with the blast-wave dynamics for transverse expansion which are incorporated 
in the Therminator model. It also indicates that a finite non-zero pion emission duration 
and outside-in burning of the source are needed to describe the data. 

5. CONCLUSIONS 

The model-independent 3D source imaging technique, involving decomposition of the 3D 
correlation function into Cartesian surface-spherical harmonics, has been applied to pion 
pairs from Pb+Pb collisions at y/s NN = 1^.3 GeV. Prominent non-Gaussian tails have been 
observed in the outward direction of the pion pair transverse momentum and in the longitudi- 
nal direction of the beam. The extracted source image is well described by the Therminator 
model which incorporates Bjorken longitudinal expansion and blast-wave transverse-flow 
dynamics. The data is consistent with a proper emission duration of 3.7 fm/c for pions at 
yfi NN = 17.3 GeV. 
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Figure 1. (color on line) Angle averaged correlation function (top panel) and source function 
(middle) and radial probability density (bottom). Filled circles show correlations from a direct 
angle-averaging of the data. Open circles represent correlations from fitting the angular 
decomposition to the data. Squares show the imaged source and correlations corresponding to the 
imaged source. The dotted and solid lines represent, respectively, the fitted Gaussian and Hump 
sources and their corresponding correlation functions (see text). 
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Figure 2. (color online) Correlation moments of multi-polarity I — 2 (left panels), and I = 4 (right 
panels) for tt + tt + and ir~ir~ pairs. Results are represented in the same way as in Fig. 1, with error 
bars indicating statistical errors for the data. Systematic errors are smaller than the data points. 
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Figure 3. (color online) Correlation C(q) (left panels) and source S(r) (right panels) function 

profiles for tt + tt + and 7r~7r~ pairs in the outward x (top panels), sideward y (middle) and 
longitudinal z (bottom) directions. The use of symbols is analogous to that in Fig. 1. Error bars 
indicate statistical errors. Systematic errors do not exceed the symbol size. Here, I = 4 moments 

make negligible contributions. 
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Figure 4. Source function profiles, S(ri), along x (a), y (b) and z (c) directions, from imaged 
correlation function (squares) and from the Therminator model (circles). 



